Appearance of Gauge Fields and Forces beyond the adiabatic approximation 
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We investigate the origin of quantum geometric phases, gauge fields and forces beyond the adia- 
batic regime. In particular, we extend the notions of geometric magnetic and electric forces discov- 
ered in studies of the Born-Oppenheimer approximation to arbitrary quantum systems described 
by matrix valued quantum Hamiltonians. The results are illustrated by several physical relevant 
examples. 
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A physical system can never be considered as com- 
pletely isolated from the rest of the universe. For a slow 
(adiabatic) cyclic variation of its environment, the wave 
function of the quantum system gets an additional geo- 
metric phase factor, known as the Berry phase [1]. In 
fact, the driving environment, the 'heavy' or 'slow' sys- 
tem, is also subject to back reaction from the 'light' or 
'fast' system. In the context of the Born-Oppenheimer 
theory of molecules, the back reaction of the light system 
leads to the appearance of a gauge field in the effective 
Hamiltonian for the slow one (the environment) ^]^] 4]. 
The gauge field consists of a vector and a scalar potential 
and turns out to depend on a quantum geometric tensor 
. It can both induce interference phenomena and mod- 
ify the dynamics through geometric Lorentz and electric 
forces 00]. 

In this note we investigate the origin of quantum gauge 
fields and forces in a more general context, by considering 
the diagonalization of an arbitrary matrix valued quan- 
tum Hamiltonian. To be precise, by diagonalization it 
is meant the derivation of an effective in-band Hamil- 
tonian made of block-diagonal energy subspaces. For 
that purpose we use the results of a powerful method 
developed recently This approach based on a new 
differential calculus on non-commutative space, where 
ft plays the role of running parameter, leads to an in- 
band energy operator that can be obtained systemati- 
cally up to arbitrary order in U. Particularly important 
for our purpose, it is possible to give an explicit effec- 
tive arbitrary diagonal Hamiltonian to order in terms 
of non-canonical coordinates and commutators between 
gauge fields (Eq. (79) in 0). We will directly apply 
this result to systems whose Hamiltonian has the sim- 
ple form iJ = T(K) -f T^(Q) and the components are 
assumed to fulfill the canonical commutation relations 
[(5i,ifj] = iMij. We then discuss how gauge fields arise 
in physical situations as various as Dirac and Bloch elec- 
trons in electric fields or Born-Oppenheimer theory. Note 
that there exists another totally different method of di- 
agonalization in a formal series expansion in h which 
uses symbols of operators via Weyl calculus To 



our knowledge this method was only applied to a Born- 
Oppenheimer- type Hamiltonian Q. 

Our approach reveals the appearance at order \}} of a 
scalar gauge potential expressed in terms of two tensors. 
One is the quantum metric tensor [5|[lol|. and the other 
is a new tensor generalizing an additional term found in 
Q for the Born-Oppenheimer case. Another very impor- 
tant consequence of the Hamiltonian diagonalization is 
the appearance of gauge invariant intraband coordinates. 
The advantage of using these coordinates is that the di- 
agonal Hamiltonian is also gauge invariant. Moreover, 
these coordinates fulfill a non-commutative algebra which 
strongly affects the dynamics through a Lorentz term and 
the gradient of a new scalar potential, generalizing thus 
the adiabatic dynamics of the Born-Oppenheimer theory. 

Hamiltonian diagonalization. Consider a Hamiltonian 
of the form 



= T(K) + y(Q) 



(1) 



where we assume that H has a matrix representation 
with T{K) non-diagonal and V{Q) diagonal. In by 
considering ft as a running parameter, we relate the in- 
band Hamiltonian UHU^ = e (X) and the unitary trans- 
forming matrix U (X) (where X = (Q, K)) to their clas- 
sical expressions through integro-differential operators, 
i.e. e (X) = 6 (eo (Xo)) and U (X.) = N {Uq (Xo)), where 
matrices with the subscribe correspond to operators re- 
placed by classical commuting variables Xo= (Qo,Kg) . 
The only requirement of the method is therefore the 
knowledge of Uq (Xo) which gives the diagonal form 
So (Xq) . Generally, these equations do not allow to find 
directly s (X), U (X), however, they allow us to produce 
the solution recursively in a series expansion in h. With 
this assumption that both e and U can be expanded in 
power series of ft, we determined in Eq. (79) of 3, the 
explicit diagonalization of an arbitrary Hamiltonian to 
order ft^. The expression of the effective n-th in-band 
energy e„ greatly simplifies for Hamiltonian given by Eq. 
([T]) . Indeed, the first order is easily obtain by an uni- 
tary transformation Uo{K.) diagonahzing T(K) giving 
UoHU+ = eo(K) + y(Q + ftA) with A = iUoVKUf^. 
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Then V is now non-diagonal. The diagonalization at 
the next order is done by an unitary transformation ma- 
trix U = Uq + HUqUi with the antihermitian matrix 
{Ui)mn = j;^Xtp|^(A)„„.Vy, which removes the off- 
diagonal elements of A and leads to corrections of order 
h"^ such that e„ in Eq. (79) of [§] becomes 

en (K, q„) = eo,n (K) + V (q„) + hH,, (2) 

where the geometric scalar potential is 

$„(Q, K) = ^d,d,V{Q)+T:,^d,V{Q)d,V{Cl) (3) 
with two gauge invariant tensors GJf and T^^ defined as 

G'^iK) = 1 5^^^^ {{A)nm{An,nn + h.C.) (4) 

and 



^'(K) = I E ^ ( + h.C.) (5) 

The tensor GJf is known as the quantum metric ten- 
sor [il[10] and T^J' is a new tensor generalizing an addi- 
tional term found in Q for the Born-Oppenheimer the- 
ory. In Eq. ([2]) we introduced the intraband coordinate 
Qn = Q + hsin where a„ is a gauge connection usually 
called the Berry connection defined as a„(K) = (A)„„ = 
i {n IVk'T')- Here \n) are the eigenstates of the non- 
diagonal part of H, i.e., T(K) \n) = eo,n(K) |n). The 
introduction of the non-canonical coordinate q„ essen- 
tial to maintain the gauge invariance of the Hamiltonian, 
implies non-canonical commutation relations [(7^,9^1] = 

z/i^ejf- with e;f (K) = |^ - |^ + [a\^,ai] the Berry 
gauge curvature in the n-th eigenstate. The Heisenberg 
equations of motion to the second order in h are 

q„ = VK£n- J(Kx0„-0„xK), K = -Vq„y (6) 

where we introduced the "magnetic field" 9*^ = e*j,9^j'^. 
The dynamics of the intraband operators leads directly 
to a Lorentz-type term. The scalar potential is a con- 
sequence of transitions between eigenstates and impacts 
the dynamics through its gradient. Working with the 
non-canonical coordinates is a short-cut to determine the 
dynamics of a system prepared in an eigenstate of the full 
Hamiltonian. This state will evolve in the same energy 
subspace as far as we can neglect higher contributions 
in the expansion in h. In comparison, the equations of 
motion derived from the Hamiltonian H do not seem to 
include a Lorentz force, and the determination of the 
" eigendynamics" can be a very difficult to achieve. An 
appealing example is given in where the "exact" 

slow motion of a massive neutral particle coupled to a 
spin is compared with the Born-Oppenheimer theory. 



We underline that our diagonalization does not need 
the adiabatic assumption, because it is an "exact" diag- 
onalization. However, the expansion in h breaks down 
in regions of mode conversion where Sq,™ — £o,n << h 
or for large values of {A^)nm — i {m \din) . Now, if in 
a particular regime the probability transition between 
different eigenstates is very small, on can neglect with 
a good approximation the off-diagonal elements. This is 
usually considered as an adiabatic approximation and we 
see that it coincides with the semiclassical approxima- 
tion. In a mode conversion region, one can easily gener- 
alize the diagonalization of to a block-diagonalization 
allowing transitions between eigenstates inside the block. 
In this case, we can therefore consider the semiclassical 
limit without having adiabaticity. 

Born-Oppenheimer approximation. Consider the fol- 
lowing Hamiltonian describing a fast system in interac- 
tion with an external environment 
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^(R,r) 



(7) 



where the fast system is described by a set of dynamical 
variables (r, p) and the slow one by coordinates (R, P). 
As in Q we consider a general kinetic energy with B, a 
positive definite inverse mass tensor. Applying the pre- 
vious results with the mapping Q — > P, K ^ R (and 
«Vk — iVr) we have V{Q) BijPiPj/2 and we ob- 
tain the following eigenvalues for the slow system (as- 
suming a non-degenerate spectrum for the fast system) 



hHn + EniK) 



(8) 



where p„ = P — ih {n |VRn) and I't.(R)) is the eigenstate 
of the fast Hamiltonian with energy i?„(R). The scalar 
potential Eq. ^ then becomes 

$„(R, P) = ^^^l^B,,+T,\'{R)BuB,kP'P'' (9) 

with the quantum metric tensor (R) — 
ReE™^„(9.n|m)(m|9,n) and y^(R) 

The term^ G^idOB,, is 
the usual part of the scalar potential discussed in 
several circumstances [llQQi whereas the term 
Ti3{Il)BiiBjkP^P^ was found in 0. Here we see 
that the Born-Oppcnhcimcr theory can be obtained 
straightforwardly from our Hamiltonian diagonalization 
to order h^. In the same manner from Eq.® we imme- 
diately get the Born-Oppenheimer equations of motion 



(R X 0, 



0„ X R) 



with Ri — Bijpl^. Similar equations of motion for 
a classical system consisting of a classical magnetic 
moment interacting with an inhomogeneous magnetic 
field 0] were studied in details. It was found that the 
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Lorentz force results from a slight misalignment of the 
magnetic moment relative to the magnetic field. This 
corresponds to the semi-classical approximation. The 
electric force is a time average of a strong oscillatory 
force induced by the precession of the magnetic moment. 
This is a kind of zitterbewegung effect. 

Particle in a linear potential. Another interesting rel- 
evant situation concerns a particle in a linear poten- 
tial exemplified here by a Bloch electron in an con- 
stant external electric field (see also ref. [ll|). Consider 
H = Ho{P, R) -I- (/?(R) with Hq the energy of a particle 
in a periodic potential and (fiCR) = — eE.R the external 
electric perturbation (and e < the charge). Using the 
mapping Q — > R, the scalar gauge potential reduces to 
<i>„(k) = e'^T^^EiEj, and the energy eigenvalues are 

e„ = eo,n (k) - eE.rn-He^^i^ry (k) E,E, (10) 

with eo,n (k) is the n-th energy band and k the 
pseudo-momentum. The intraband position operator 
is r„ = R + ih{un iVkUn) with |u„(k)) the peri- 
odic part of the Bloch wave function and T'^-'(k) = 
Re J^m^jin ^^'""^"'"^ig" _ Introducing the " magnetic 

■E.Vk^'n the 



field" w,, 



eE' 



^ ExVk$„ and x„ 



equations of motion are 

h 



eE' 



(k X f]„ -I- ri„ X k) + h^Xn^^ k = eE 



where r2„ = 0„ -I- ?iu;„ . This shows that $„ contributes 

to the Lorentz term ftk x ri„ known as the anomalous 
velocity which is orthogonal to the applied electric field. 
This anomalous velocity is at the center of many recent 
experimental and theoretical works. <&„ contributes also 
to the velocity in the direction of E, through the term 

Berry phase. The linear potential case has another 
interest. It allows us to also consider the fast system 
and derive the Berry phase in a different way. Indeed, 
consider a time dependent Hamiltonian H(t) and intro- 
duce the differential operator D = H (t) — Pq where 
Po = ihd/dt is the conjugate of time which is treated 
formally as an operator such that [Pq: t] = The time 
dependence is due to the time evolution of some param- 
eters x(t) describing the environment. To transform the 
system of differential equations (Schrodinger equation) 
D \'${t)) = 0, which couples all components of |^'(i)) into 
a decoupled set of differential equations, we introduce an 
unitary transformation (t)) = U{t) |^ (i)) such that 
U{t)D{t,Po)U+{t) = A(i,Po) is a diagonal differential 
operator and A{t, Pq) {t)) = 0. Therefore the time 
evolution is given by 1^-' (t)) = (o)) . Since 

A(t) = K{t) + Pq is diagonal, no time ordered product is 
required. Returning back to the initial state we have 

1^- [t)) = U+{t)e^ ■I'o ^<-*^<^'U{0) I* (0)) (11) 



A system prepared in a state |A„(0)) which is an eigen- 
state of D, i.e., D(0) |A„(0)) = A„(0) |A„(0)) will evolve 
with A(t) and thus stays in the instantaneous eigenstates 
\An{t)) of D{t) (for simplicity we assume non degen- 
erate eigenvalues). In this case the wave function be- 
comes l'^ (t)) = e^-/^o''^"(*)'^* |A„(i)). Since eigenstates 
of D instead of H are considered, the time evolution Eq. 
()11|) is non-adiabatic. In general we need an approxi- 
mation scheme for the diagonalization of D and we will 
use the expansion to order fi^. The problem of finding 
An = An (t) -t- Po is formally equivalent to the Bloch elec- 
tron example discussed above with K t, Q = R ^ Pq 
and eE ~ 1. We obtain from Eg. flOl) 



A. 



en{t)-ih{n n) +h'^Tn (t) 



(12) 



where £„ (t) and n{t) are instantaneous eigenvalues and 
eigenstates of H, i.e., H{t)\n{t)) — e„(t)|n(i)), and 

(n\m){m\n\ 

Tn (t) = ^^zZm^n e -e ' Therefore for a periodic 
motion of period T , not necessarily adiabatic, such that 
|A„(T)) — |A„(0)) (single valued eigenstates), we have 



with 



|*(T)) 



£n{t)dt -f i 





l*(0)) 








71^ dt +h J 







(13) 



Tn {t) dt 



The phase 7„ appears as an expansion in power of h. 
The first term is the usual dynamical phase and the sec- 
ond one the geometric Berry phase independent of h and 
of the velocity of parameters x{t). The additional phase 
h /g Tn (t) dt of order H is non-geometric as it depends 

on X. It cancels in the infinitely slow a; ^ adiabatic 
regime, which thus coincides with the semiclassical ap- 
proximation. Diagonalization at order thus goes be- 
yond the adiabatic approximation and takes into account 
transitions between eigenstates. Quantitatively, if the 
system is prepared in an eigenstate |n(0)) of H{Q), then 
\'^{t)) is given by Eq. with U ^ Uq + KUiUq where 



Ui {t\ 



■ - e^it)-eAt) 

expansion up to order h: 



, so that we have the following 



|4-(t)) = i'o'^-W''* |n(t)) + ^ {e^^°^-^'^'^'A„,n{t) 



The magnitude of transitions is then controlled by the 

i/m(t)|n(t)\ 

term Amn = g {t)-e (t) '^hich is neglected in the adi- 
abatic limit (rn{t) n{t)^ 0. Note that [*(<)) is nor- 
malized to unity at order h only. A normalization at a 
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higher order needs an expansion of U to the same order 
i- 

In principle deviation from adiabaticity could be mea- 
sured by interferometry. Consider a periodic two states 
system, and write the initial state in the eigenbase 
|n(0)) = |A„(0)) + M„,„(0) |A™(0)) . Then, after one 
cycle l^- (T)) = e"^^- |A„(0)} + fe-'T"A„™(T) |A,„(0)) . 
For an observable O which does not commute with H 
one will find in the average ('3/ (T)| O (T)) an interfer- 
ence term 2hRe (A„„(r) (A„(0)| O |A„(0)) e"^*'^"-'^'")) 
which would signal deviation from adiabaticity. As dis- 
cussed below, this interference effect is formally equiva- 
lent to the zitterbewegung of Dirac particles. 

Dirac particle in an external potential. 

We will now show that our formalism can also be used 
for relativistic Dirac particles, which are usually treated 
with the Foldy Wouthuysen approach [l^]. The Hamil- 
tonian is (with c = 1) 

H = a.p + (3m. + V (R) (14) 

where a. and /3 are the usual (4 x 4) Dirac matrices and 
V (R) is the external potential. The matrix diagonalizing 
the free part of the Hamiltonian ?7o (a.p + (im) Uq = 
(3E with E = -I- is the usual Foldy Wouthuysen 
unitary transformation Un — , '^'tl?^^^f/2 • For the Dirac 
particles we have two energy subspaces £± of dimension 2 
corresponding to the positive and negative energy. Now 
with the correspondence Q ^ R, K — > p and formula 
Eq. ^ , one easily sees that the diagonal matrix can be 
written 

e (p, r) = /3£; (p) -f F (r) -f /3$ (15) 

The position operator is given by the (4 x 4) matrix 
r — R+ 2g(^-Hn.) "^ith S = 1 ® (T where cr are the Pauli 
matrices. The band index of the scalar potential has been 
transferred to the matrix /3, and we have G*^ = j^g^-' 

and T'^ = g^.gy with the notation g'^ = 5'^ - so 
that finally we can write 

'^-^9''{^^^,V+^^,V^,V) (16) 

If for central potential one can neglect the contribution 
-^diVdjV, this is not always true and for some potentials 
both terms in Eq. can be of the same magnitude. In 
fact for constant electric field V — — e£.R, the first term 

2 1-2 

vanishes and <I> = ^^g^^ £i£j. 

In the non-relativistic limit p << m, $ becomes 
$ « {/S.V+^{VVf)) + 0{h^p^/m'^) which gives 
two contributions. The first one is the usual Darwin 

i-2 

term ^^AV traditionally obtained as the result of the 
Foldy Wouthuysen transformation expanded in power of 
1 /m. The second term C^V) of higher order in 1/m 



is usually not considered in the Foldy Wouthuysen ap- 
proach. It is also interesting to note that the external 
potential in the non relativistic limit can be expanded as 
F (r) « y (R) + j^-E. {VV X p) + 0(SV/™'') where 
j^'E- (VV X p) is the spin-orbit coupling term. There- 
fore the Hamiltonian can be approximated as 

/ p2 p4 \ 

A Born-Oppenheimer treatment of the Dirac equation 
where the spin is the fast variable and the momentum 
the slow one has led to the same Hamiltonian Eq. pT]) 
but without the scalar potential [l^ . This corresponds to 
the scmiclassical approximation. The additional electric- 
type potential $ is a consequence of transitions between 
energy levels. This is in agreement with the usual in- 
terpretation of the physical origin of the Darwin term, 
the zitterbewegung phenomenon, whereby the electron 
does not move smoothly but instead undergoes extremely 
rapid small-scale fluctuations due to an interference be- 
tween positive and negative energy states. 
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